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El resumen sera traducido al espanol por los editores. We summarize results on the thermal and 
• ionization structure of self-similar, magnetically-driven, atomic disk winds heated by ambipolar diffusion. We 
improve upon earlier work by Safier by considering (1) new MHD solutions consistent with underlying cold 
keplerian disk equilibrium, (2) a more accurate treatment of the micro-physics, and (3) predictions for spatially 
resolved forbidden line emission (maps, long-slit spectra). The temperature plateau ~ 10 4 K found earlier 
is recovered, but ionization fractions are revised downward by a factor of 10, due to previous omission of 
Oh' thermal speeds in ion-neutral momentum-exchange rates. The physical origin of the temperature plateau is 
outlined. Predictions are then compared with T Tauri star observations, with emphasis on the necessity of 
suitable beam convolution. Jet widths and variations in line profiles with distance and line tracer are well 
reproduced. However, predicted maximum velocities are too high, total densities too low, and the low-velocity 
[O i] component is too weak. Denser, slower MHD winds from warm disks might resolve these discrepancies. 



•rH . ABSTRACT 

We summarize results on the thermal and ionization structure of self-similar, magnetically-driven, atomic disk 
winds heated by ambipolar diffusion. We improve upon earlier work by Safier by considering (1) new MHD 
solutions consistent with underlying cold keplerian disk equilibrium, (2) a more accurate treatment of the 
micro-physics, and (3) predictions for spatially resolved forbidden line emission (maps, long-slit spectra). The 
temperature plateau ~ 10 4 K found earlier is recovered, but ionization fractions are revised downward by a 
factor of 10, due to previous omission of thermal speeds in ion-neutral momentum-exchange rates. The physical 
origin of the temperature plateau is outlined. Predictions are then compared with T Tauri star observations, 
with emphasis on the necessity of suitable beam convolution. Jet widths and variations in line profiles with 
distance and line tracer are well reproduced. However, predicted maximum velocities are too high, total 
densities too low, and the low-velocity [O i] component is too weak. Denser, slower MHD winds from warm 
disks might resolve these discrepancies. 

Key Words: MAGNETOHYDRODYNAMICS — ISM: JETS AND OUTFLOWS — STARS: PRE-MAIN- 
SEQUENCE 

1. INTRODUCTION 
Collimated mass ejection in young T Tauri stars 
(TTS) is intimately correlated with the accretion 
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process (Cohen et al. 1989; Cabrit et al. 1990; Harti- 
ganetal. 1995; hereafter HEG95). It is currently be- 
lieved that magnetic forces are responsible for both 
the high ejection efficiency (Mj/M acc ~ 0.01 — 0.1; 
Hartigan et al. 1994, HEG95) and the high degree 
of collimation of these winds. Yet, the exact field 
structure and flow dynamics arc not well established. 
As a result, the role played by magnetized ejection 
in the physics of TTS, in particular in extracting an- 
gular momentum from the circumstcllar disk and / or 
stellar surface, remains a major enigma. 

Ejection from a wide range of disk radii has been 
invoked to explain the compact low-velocity [O i] 
emission peak observed in accreting TTS (Kwan & 
Tademaru 1988, 1995; Hirth et al. 1994; HEG95; 
Kwan 1997). However, it is still unclear whether the 
high-velocity "microjets" observed towards several 
TTS could trace the outer collimated regions of the 
same disk wind, or whether they require a distinct 
ejection component. 

A new tool to discriminate between models con- 
sists in confronting theoretical predictions with re- 
cent spatially resolved observations of the inner wind 
structure of TTS in optically thin forbidden lines of 
[O i], [S n], and [N n] (see Bacciotti et al., Dougados 
et al., Ray et al., this volume, and refs. therein). 

Of the wide variety of models available in the lit- 
erature, only two classes have self-consistent station- 
ary MHD solutions that have been used for detailed 
observational predictions: (1) 'disk winds' from a 
wide range of disk radii (Blandford & Payne 1982; 
Ferreira 1997) and (2) 'X-winds ' from a tiny region of 
open stellar field lines near the disk corotation radius 
(Shu et al. 1995). Both are in the regime of "cold" 
MHD, where ejection is only possible by magneto- 
centrifugal launching on field lines sufficiently in- 
clined from the disk axis, as thermal energy is in- 
sufficient to cross the potential barrier (unlike in the 
solar wind). 

Synthetic maps and long-slit spectra were first 
presented by Shang et al. (1998) and Cabrit et al. 
(1999) for X-winds and disk winds respectively, using 
a parametrized temperature and ionization fraction. 
More reliable predictions require actually solving for 
the wind thermal and ionization state, given some lo- 
cal heating mechanism. Shock heating was recently 
invoked to explain line ratios in the outer regions of 
microjets (Dougados et al., Lavalley-Fouquct et al. 
2000) but this mechanism is not fully consistent with 
a stationary solution and involves additional free pa- 
rameters. Other modes of mechanical energy dissi- 
pation (e.g. turbulence) have the same drawbacks. 

Here, we consider only a heating process intrinsic 



to MHD winds of low ionization, namely ion-neutral 
frictional drag ( 'ambipolar diffusion heating ') , which 
was first applied to predict integrated forbidden line 
emission in disk winds by Safier (1993a, b). It re- 
quires no extra free parameter and is self-consistently 
determined by the disk wind MHD structure. It will 
yield the "minimum" possible emission fluxes and 
ionization level, but emission maps and line profiles 
should still provide useful tests of the model collima- 
tion and dynamics. A detailed account of our results 
is given in Garcia et al. (2001). 

2. THERMAL AND IONIZATION STRUCTURE 
2.1. Method and improvements over previous work 
We follow the same general approach as in the 
pioneering work by Safier (1993a). We follow a fluid 
cell along a given flow streamline, characterized by 
zuq, the field line footpoint radius in the disk mid- 
plane, and solve for its coupled temperature and ion- 
ization evolution, including all relevant heating and 
cooling terms. We check a posteriori that the sin- 
gle, cold fluid assumptions made in computing the 
dynamical MHD solution remain valid, i.e. that 
drift speeds between ions and neutrals are small 
compared to the bulk flow speed, and that ther- 
mal pressure gradients are negligible compared to 
other forces. Our equations, heating/cooling terms, 
and self-consistency checks are described in Garcia 
et al. (2001). Important improvements to the work 
of Safier (1993a) include: 

- Use of new self-similar cold disk winds solu- 
tions from Ferreira (1997), where vertical equilib- 
rium of the keplerian accretion disk is treated self- 
consistently. The mass loading and magnetic lever 
arm are then fixed by a single parameter: the ejec- 
tion efficiency in the disk £ = d In M acc /d In vjq. 

- Treatment of the ionization evolution of heavy 
elements (C, N, O, S, Ca, Mg, Fe, ...), photo- 
ionization heating, and cooling by hydrogen recombi- 
nation lines, using the MAPPINGS Ic code (Binette 
et al., 1985, Binette & Robinson 1987). 

- Correct computation of ion-neutral momentum 
exchange rates for ambipolar diffusion, including 
thermal speeds (cf. Draine 1980). 

- Depletion of heavy atoms into grains outside the 
dust sublimation cavity (Savage & Sembach 1996). 

2.2. MHD solution parameters 
To limit the number of free parameters, we fix 
the disk aspect ratio s = h/r = 0.1, as indicated by 
observations of HH 30 (Burrows et al. 1996), and 
set the magnetic turbulence level parameter (con- 
trolling disk accretion across field lines) to a m = 1 
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Fig. 1. Temperature T, electronic density n e , and proton 
fraction f p — n(H + )/n,H along two streamlines anchored 
at roo = 0.1 AU (left column) and wo = 1 AU (right 
column), for a disk wind with £ = 0.007. For compar- 
ison, dashed curves show calculations with local ioniza- 
tion equilibrium imposed. M acc increases in the direction 
of the arrow from 10 -8 and 10 -5 M@ yr _1 by factors of 
10. Dotted vertical lines mark the locus of the Alfven 
surface (x — 10) and the point of maximum streamline 
opening (x — 2 x 10 4 ) before refocusing 



(cf. Ferreira 1997). The central star mass is taken 
as 0.5 Mq. Thus only 2 free global parameters are 
left: 

- The ejection efficiency parameter £, which typ- 
ically ranges between 0.005 and 0.01 for the above e 
and a m . Our disk wind extends from wi to 10 x Wi 
and thus £ ~ Mj Ct /M acc as observed in TTS. 

- The disk accretion rate, M acc , which deter- 
mines both the jet density scalings and the pho- 
toionizing UV flux from the accretion disk boundary 
layer. We consider values between 10 -8 and 10~ 5 
Mq yr _1 , covering the observed range in accreting 
TTS (HEG95). 




10 1 10 10 10 10 
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Fig. 2. Heating and cooling terms (in erg s 1 cm 3 ) for 
the same model as in Fig. 1. Top: M acc =10 _6 MQyr _1 , 
Bottom: M acc =lO~ 7 A/ yr _1 . 



2.3. Results 

Figure [l] presents the computed temperature, 
electronic density n c , and proton fraction / p along 
two streamlines for a typical solution with an inter- 
mediate £ value of 0.007 and various mass accretion 
rates. The main results are the following: 

- Temperature reaches a plateau around 10 K 
over most of the jet extent for the range of parame- 
ters applicable to T Tauri disk winds. 

- The proton fraction / p rises steeply at the 
wind base, then freezes out in the far jet region 
(x = z/zuq > 300) where densities are low and ion- 
ization timescales exceed dynamical timescales. In 
the temperature plateau, / p is roughly inversely pro- 
portional to Macc an d n7o. 

- The electron density n c is dominated by pho- 
toionization of heavy elements at the wind base, but 
is weakly dependent on M acc further out. 

- The main heating process is ambipolar diffu- 
sion, Fjjrag, and the main cooling term is adiabatic 
expansion, A ac ii a - A close match between the two 
terms is established in the temperature plateau re- 



gion (see Fig. 2.3). 

Remarkably, the same behaviors were found for 
the solutions investigated by Safier (1993a). In the 
next section we explain why they are generic prop- 
erties of self-similar cold MHD disk winds. 
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Fig. 3. Left: F(T) in erg g cm s _1 assuming local ionization equilibrium and photoionization of Na only (solid) or 
C (dashed). Center: G(x) in erg g cm 3 s" 1 for models A, B, C (bottom to top), uj = 0.1 AU and M acc = 1CT 6 M© 
yr _1 . Right: Temperature for model B assuming ionization equilibrium (dashed), compared with T c (x) as given by 
Eq. fl] (dot-dashed). zu = 0.1 AU and Af acc = 10" 8 to 10~ 5 M©yr _1 , from top to bottom. 



2.4. Physical origin of the hot temperature plateau 

It is useful to write the temperature evolution 
equation as a function of \ — z/wq as: 



d\nT 
d\nx 



2 d In n 



3 d\nx 

s- 1 



x (G(x) 



-A-adia 

- 1 



(i) 



where h is the wind particle density, 6 is a positive 
function of order 1 before jet recollimation, 



G(x) 



fi 2 (v ■ V)n 



(2) 



is another positive function (before recollimation) 
that depends only on the MHD wind solution, while 

F(T) = kT(l + / e )m in /i/ n < a m v > x (-B-Y (3) 

depends only on the local temperature and ionization 
state of the gas. 

One can see that F(T) always increases as a func- 
tion of temperature: either linearly oc T /; when ion- 
ization is fixed (by photoionization or freeze-out), or 
much more steeply when collisional ionization of H is 
efficient, around 10 4 K. These two regimes of F(T) 
are illustrated in Fig. || in the case of local ioniza- 
tion equilibrium. The function G{\) is a l so plotted 
in Fig. ^ for our solutions. It increases very rapidly 
at the base of the flow and then stabilizes in a plateau 
beyond the Alfven point. 

Now let us define a temperature T c (x) such that: 



F(T C ) = G( X ) 



drag 



= A 



adia- 



(4) 



Since F increases with temperature, it can be read- 
ily seen that if at a given point T > T c (x), then 
G(x)/ F(T) < 1 and the gas will cool. Conversely, if 
T < T c (x), the gas will heat up. Thus, the fact that 
F(T) is a rising function introduces a feedback that 
tries to bring the temperature near its local equilib- 
rium value T c (x), and A adia near r drag . 

However, T ~ T c (x) is not necessarily a possible 
solution of Eq. 0: in that case the right-hand term 
is close to zero (by definition of T c in Equ. |]) hence 
one must also have d\nT c / d\nx <C 1 for consistency: 
Only when T c (x) is flat can T converge to T c , and can 
we have a temperature plateau with A a d; a — rd rag . 

Indeed, T c is, for our models, a flat function. We 
can understand this by differentiating Eq. 0: 



d\uG 
G?lnx 



d\nF 
d\nT 



T=T C 



d In T c 
dlnx 



(5) 



Hence |dlnT c /c?lnx| <C 1 is equivalent to 
\dlnG/dlnx\ < \dln.F/dlnT\. This is fulfilled for 
our models: Below the Alfven surface, G varies 
a lot, but collisional H ionization is sufficiently 
close to ionization equilibrium that F(T) still rises 
steeply around 10 4 K (cf. Fig. ||). For our nu- 
merical values of G, we have T c ~ 10 4 K and thus 
\d\nG/dlnx\ < \dlnF/dlrxT\. Further out, where 
ionization is frozen out, we have dlnF/dhiT = 1 
(because F oc T/;) but it turns out that in this re- 
gion G is a slowly varying function of x, and thus we 
still have |dlnG/dlnx| < |dlnF/dlnT|. 

The behavior of f p when G = F can also be un- 
derstood: in self-similar disk wind models, G scales 
as l/(n7 Af acc ). Since F(T) ex Tf p (when / e < 1) 
and / p is a high power of T close to ionization equilib- 
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rium, most of this scaling is absorbed by / p , while T 
remains around 10 4 K for a wide range of model pa- 
rameters. Note also that inclusion of thermal speeds 
(cf. Draine 1980) increases < ai n v > by ~ 10 com- 
pared to Safier (1993a) and decreases our values of 
f p by the same amount, since only the product of 
these two terms enters in F(T). 

We conclude that MHD winds heated by ambipo- 
lar diffusion have a hot temperature plateau only 
when several conditions are met: (1) G must be such 
that F = G happens around 10 4 K, (2) the wind 
must be in ionization equilibrium or near it in re- 
gions where G is a fast function of \\ (3) once we 
have ionization freezing, G must vary slowly. For 
example, in the models of Ruden et al. (1990), ion- 
ization is quickly frozen while G oc 1/r, hence no 
temperature plateau can be established (one has in- 
stead T oc 1/r). 

3. PREDICTED FORBIDDEN LINE EMISSION 
AND COMPARISON WITH OBSERVATIONS 

We compute observational predictions assuming 
that the disk wind extends from 0.07 AU (typical 
disk corotation radius for a T Tauri star) to 1 AU. 
Inside corotation, the disk is expected to be trun- 
cated by the stellar magnetosphere, while beyond 
1 AU, the wind should become molecular and not 
contribute to forbidden line emission (Safier 1993a). 

Given the strong gradients in physical conditions 
present in the wind, care was taken to convolve syn- 
thetic maps and spectra to typical resolutions, to al- 
low meaningful comparison with observations. The 
successes and failures of our model predictions are 
outlined below. 

3.1. Wind morphology: Emission maps 

The predicted maps for our models (Fig. ^) suc- 
cessfully reproduce two features present in observed 
microjets: (1) an unresolved peak slightly shifted 
from the stellar position, (2) jet-like emission appear- 
ing collimated within 200 AU of the star. Note that 
the MHD disk wind solutions considered by Safier 
1993a did not collimate over such scales. 

The shift of the unresolved peak cannot be used 
as model test because observed shifts are severely bi- 
ased by accuracy of continuum subtraction (2% er- 
rors in subtraction imply ~ 4 AU error for bright jets 
(Garcia et al. 1999)) and by possible obscuration 
by a flared disk. 

The jet FWHM, in contrast, provides an excellent 
quantitative test of the model: It is almost indepen- 
dent of the inclination angle, accretion rate and line, 
but is strongly dependent on the MHD dynamical 
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Fig. 5. Predicted jet transverse FWHM (obtained 
with gaussian fits) for models A, B, C (with £ = 
0.01,0.007,0.005) and two different beam sizes: 28 
AU and 2.8 AU corresponding in Taurus to 0.2"and 
0.02" (VLT diffraction limit). Note the strong depen- 
dence on £ and the bias introduced by the beam size. 



solution (Figs. [|]|). Note that it is also severely bi- 
ased by the beam size (Fig. ||), an illustration of the 
importance of proper convolution before comparing 
to observations. 

Predicted jet FWHM for models A and B (f = 
0.007 — 0.01) agree very well with observations of 
T Tauri microjets at the same resolution (Douga- 
dos et al., 2000 and in this volume). Models with 
lower £ (e.g. model C) are excluded. Note that at 
sufficiently high spatial resolution, disk jets should 
appear hollow, hence it would be important to get 
very high angular resolution data. 

3.2. Wind kinematics: Line profiles 
Typical long-slit spectra are presented in the 
right panels of Fig. [|, and integrated spectra for a 
wider range of parameters are plotted in Fig. |^. The 
predictions are successful in several respects: 

- A compact low-velocity component (LVC) is 
produced near the star (Fig. ^), originating at the 
base of the slowly rotating outer disk streamlines. 
The LVC is stronger in [O i] than in [N n], and is 
stronger for low accretion rates (Fig. as observed 
(HEG95). 

- An extended high-velocity component (HVC) 
is also present, tracing the accelerated regions of the 
wind. It is displaced further from the stellar position 
than the LVC, and the displacement is larger in [N n] 
than in [O i] (Fig. g), as observed (Hirth et al. 1997). 
The [N n] line profile peaks at the blue edge of the 
[O i] line profile (Fig. ||), as also observed (HEG95). 

However, we fail to reproduce some observed fea- 
tures: 
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Synthetic [sii] maps Long-slit spectra along jet axis 

Model A Model B Model C 

^=0.01 ^=0.007 ^=0.005 Model A ^=0.01 Model B ^=0.007 Model C ^=0.005 
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Fig. 4. Emission maps (left) and long slit spectra along the jet axis (right) for the all models and several forbidden 
lines. The inclination angle is i = 60° from pole-on, and mass accretion rate is M acc = lO -6 Af0yr -1 . Intensity maps 
are convolved with a 28 AU resolution beam (0.2"in Taurus) and long-slit spectra are convolved by 70 AU x lOkm/s, 
representative of current ground-based spectroimaging performances. Contours decrease by factors of 2. 



- The relative intensity of the LVC with respect 
to the HVC is too small. At i < 60°, it does not 
appear anymore as a distinct component in the [O i] 
profile. This weakness of the LVC stems from the low 
temperature and ionization below the Alfven surface, 
where low-velocity gas is located. 

- Predicted centroid and maximum radial veloc- 
ities of [O i] profiles agree well with observations of 
the DG Tau microjet, but are too high compared 
with typical TTS profiles (HEG95) unless most stars 
are observed at i > 80° (see Fig. ||). Terminal veloci- 
ties in cold disk wind models are ~ Vkep^oVV?? so 
even our highest £ = 0.01 (model A) gives excessive 
speeds. 

3.3. Integrated line fluxes 

Several trends are well reproduced: Predicted in- 
tegrated fluxes increase linearly with accretion rate, 
as observed (HEG95). This is an interesting con- 
sequence of the weak dependence of n c in the jet 
on M acc (cf. Fig. 1): emissivity (oc n c nu) is then 
proportional to total density njj, instead of n 2 H as 
would be the case for a constant ionization fraction. 
We also find a correlation between integrated fluxes 
in [O i] and [S n] , with a slope close to that observed 
(see Cabrit et al. 1990). 

Quantitatively, however, integrated fluxes are 
systematically too weak compared to observations 
of TTS, for the same range of accretion rates, i.e. 
10~ 8 - 10~ 5 M Q yr- 1 (HEG95). We find L([0 i]) in 
Lq ~ M acc in Mq yr _1 . The discrepancy is typically 
a factor 30 for a given M acc . In the next section, we 



show that the flux deficit stems from insufficient n e , 
but also possibly from insufficient total density nn 
in our models. 

3.4. Forbidden line ratios and total wind densities 

As noted by Bacciotti & Eisloeffel (1999), forbid- 
den line ratios of [S n] 6716/6731, [N n]/[0 i], and 
[S n]/[0 i] reflect directly the values of n e , f c , and 
T. Hence, they test mainly the heating mechanism, 
not the underlying dynamical solution. 

Spatially resolved line ratios within 200 AU of 
the star are available for four TTS microjets: HH 30, 
Th 28, DG Tau, and RW Aur (Bacciotti & Eisloeffel 
1999; Lavalley-Fouquet et al. 2000; Bacciotti et al. 
and Dougados et al., this volume). Comparison with 
our predictions, at the same beam resolution and dis- 
tance from the star, indicates that our models have 
slightly too high temperatures, and beam-averaged 
n e and f e that are too low by up to a factor 30. Two 
important remarks are in order: 

(1) because of volume effects and strong density 
gradients, beam-averaged line ratios differ greatly 
from local line ratios, and give very little weight to 
the inner densest parts of the wind. Hence it is es- 
sential to apply proper convolution to perform mean- 
ingful comparison with observed ratios. 

(2) No ad-hoc heating rate tuning is present in 
our model. Predicted line ratios would improve if 
we include extra mechanical energy deposition, e.g. 
in shocks, as illustrated in Dougados et al. (this 
volume), or with a dissipation length prescription, 
as done in Shang et al. (this volume). 
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Changing MHD solution 
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' I 

I ' I M I K'' 


^^1^ i' 











-300 
V (km/s) 



-600 -300 
V (km/s) 



Fig. 6. Influence of model parameters on line profiles. 
Left: effect of accretion rate. The arrow points in the 
direction of increasing M acc . Middle: effect of MHD 
solution. Right: effect of inclination angle. The arrow 
points in the direction of increasing i. 
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For the latter reason, the total density n# = 
n e /f e is a more useful quantity that can be directly 
compared with model predictions independently of 
the ionization process (although it can still be af- 
fected by shock compression and by beam averag- 
ing). For accretion rates of 1CP 6 M Q yr -1 , total 
densities on the innermost streamline for our models 
are too low by a factor 10-100 compared to obser- 
vations of the above four bright T Tauri microjets. 
Shock compression could explain a factor of 10 (Har- 
tigan et al. 1994), but it seems that our models have 
intrinsically too low density close to the star. 

4. CONCLUDING REMARKS 
We extended the original work of Safier 
(1993a, b), and compute the thermal structure of a 
disk wind jet, assuming current dissipation in ion- 
neutral collisions ambipolar diffusion, as the major 
heating source. Improvements include: a) detailed 
dynamical models where the disk is self-consistently 
taken into account; b) ionisation evolution for all rel- 
evant "heavy atoms" ; c) radiation cooling by hydro- 
gen lines, recombination and photoionisation heat- 
ing; d) correct H — H + momentum exchange rates; 
and e) more detailed dust description. We still ob- 
tain warm jets 10 4 K but with ionisations fractions 
10 to 100 times smaller than Safier, due to larger 
H — H + momentum exchange rates (including the 
dominant thermal velocity contribution) and to dif- 
ferent MHD wind dynamics. 

We have presented a complete set of predictions 
for forbidden line emission, where we have stressed 
the crucial effect of beam convolution. The model 



Hot, dense jet 



Cold, tenuous jet 





Fig. 7. Upper 4 panels: New MHD disk wind solu- 
tions with large heating at the disk surface (solid curves) 
have: more collimated streamlines (top left), lower termi- 
nal velocity (top right), and higher density (middle left) 
than a model with small heating (dashed curves). Bot- 
tom panels: Local energy dissipation rate (Joule and 
viscous; r e fi), and prescribed entropy injection Q for the 
2 models displayed above. See Casse & Ferreira (2000) 
for details. 



reproduces several observed trends: (1) images show 
an unresolved peak and an extended high-velocity 
jet, of width compatible with observations. (2) Line 
profiles present a low- velocity component (LVC), 
compact and near the star, and a high- velocity com- 
ponent (HVC), tracing the jet. The systematic dif- 
ferences between [O i] and [N n] profiles are repro- 
duced. (3) Line fluxes are proportional to M acc , and 
the [O i]- [S n] correlation slope is recovered. 

Other line ratios ([S n]6716,6731, [N n]/[0 i]) are 
not well reproduced, but we stress that they trace 
only the excitation conditions, which depend mainly 
on the heating mechanism, and do not test the un- 
derlying dynamical model; a better fit could be read- 
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ily obtained with an extra tunable heating rate, such 
as done e.g. by Shang et al. (this volume) or with 
shocks (Dougados et al., this volume). 

Two intrinsic properties of our dynamical models 
are not fully successful: Wind terminal velocities ap- 
pear too high, and densities at the wind base appear 
too small, when compared with current estimates in 
bright microjets. One possible improvement would 
be to relax the assumption of a cold disk wind, and 
include heating at the disk surface (Casse & Fer- 
reira 2000). MHD disk winds with higher density 
and lower terminal velocities can then be obtained 
(see Fig. Q) . The self-consistent thermal structure of 
these solutions remains to be investigated. 
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